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ABSTRACT 

We present a mass estimate of the Planck-discovered cluster PLCK G100.2-30.4, derived from a weak lensing analysis of deep 
SUBARU griz images. We perform a careful selection of the background galaxies using the multi-band imaging data, and undertake 
the weak lensing analysis on the deep (Ihr) r-band image. The shape measurement is based on the KSB algorithm; we adopt the 
PSFex software to model the Point Spread Function (PSF) across the field and correct for this in the shape measurement. The weak 
lensing analysis is validated through extensive image simulations. We compare the resulting weak lensing mass profile and total 
mass estimate to those obtained from our re-analysis of XMM-Newton observations, derived under the hypothesis of hydrostatic 
equilibrium. The total integrated mass profiles are in remarkably good agreement, agreeing within Icr across their common radial 
range. A mass M 500 ~ 7 x 10^"^ M© is derived for the cluster from our weak lensing analysis. Comparing this value to that obtained 
from our reanalysis of XMM-Newton data, we obtain a bias factor of (1-^) = 0.8 + 0.1. This is compatible within Icr with the 
value of (1-b) obtained by [Planck Collaboration XXrV| from their calibration of the bias factor using newly-available weak lensing 
reconstructed masses. 

Key words. Galaxies: clusters: individual: PLCKGlOO.2-30.4 - Gravitational lensing: weak - X-rays: galaxies: clusters - Cosmol¬ 
ogy: dark matter 


1. Introduction 


Clusters of galaxies are a crucial probe for a variety of key cos¬ 
mological issues. Among their many applications, knowledge of 
the mass function and its redshift evolution allows constraints 
to be put on the matter power spectrum and limits to be set on 
several cosmological parameters. The construction of the cluster 
mass function generally involves the following steps: 1) clusters 
should be identified in large surveys; 2) mass proxies should be 
used to estimate their masses; and 3) the mass function should be 
constructed by counting clusters in given mass and redshift bins. 
This process relies on the fundamental assumption that the scal¬ 
ing relations between the mass and its proxies are well-calibrated 
and unaffected by selection biases (or that the biases have been 
accounted for). It has been successfully implemented both using 
X-rays surveys as the ROSAT All-Sky Survey, where proxies 
like the gas mass, the temperature, or the X-ray luminosity have 
been used to estimate the masses (e.g., jVikhlinin et al.|[2009 


Mantz et^|2010| ), and in optical photometric optical surveys 
(e.g. the SDSS), where the masses are usually derived from the 


* Based on data collected at Subaru Telescope (University of Tokyo). 


cluster richness (e.g., |Rozo et al.|2009] ). More recently, searches 
for galaxy clusters have been started via the detection of the 
Sunyaev-Zel’dovich (SZ) effect on the photons of the Cosmic- 
Microwave-Background. The newly-established millimetric ob¬ 
servatories, both ground- (South Pole Telescope and Atacama 
Cosmology Telescope) and space-based (Planck), have provided 
catalogues of hundreds to thousands of clusters. The first clus¬ 
ter detections from the Planck all-sky survey were released in 
the Planck Early SZ sample ( [Planck Collaboration VIII||2Q11| ). 
A crucial part of the validation process involved use of snapshot 
XMM-Newton observations; in the first two campaigns, 21 out 
of 25 newly identified cluster candidates were confirmed to be 
real ( [Planck Collaboration IX|20lT] ). 


The XMM-Newton snapshot observations appeared to reveal 
a class of massive clusters with low X-ray luminosity and dis¬ 
turbed morphology, which may indicate a population of massive 
clusters that are under-represented in X-ray surveys. Since the 
cluster masses in SZ surveys are estimated from the integrated 
Compton Y parameter via scaling relations calibrated with op¬ 
tical or X-ray observations, this may have an effect on the esti¬ 
mate of the mass function, and thus our understanding of cluster 
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Fig. 1. Composite 2 ri image of the inner 400" X 400" region of the 
PLCKGIOO field. 
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Fig. 2. Zoom of the inner 100" x 100" region of the PLCKGIOO field. 
Labels mark the 10 brightest galaxies. Several strong tensing gravita¬ 
tional arcs are visible. 


growth and evolution. It is therefore of fundamental importance 
to study these systems in more detail. 

Lensing is a powerful method to recover the surface den¬ 
sity field of clusters, thus allowing to trace the spatial distribu¬ 
tion of dark matter. Simulations have shown that provided wide- 
field observations of excellent quality are available, lensing al- 
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lows an accurate measu rement of the mass profile of clusters 
( [Meneghetti et al.||2Q10|), and hen ce the radial mass/luminosity 
(M/L) profile ( [Medezinski et al.||2010] ), and its dependence on 
early- / late-type galaxies in the cluster. Combined with X-ray 
observations, lensing can help to characterize the dynamical 
state of the clusters, and probe deviations from hydrostatic equi- 
librium which ma y bias the X-ray masses ( [Mahdavi et al.|200^ 
[Zhang et al.|2Q10| ). 

In this paper we describe the first weak lensing analysis of 
the cluster PLCK GlOO.2-30.4 (hereafter PLCKGIOO), based on 
deep SUBARU SuprimeCam images. PLCKGIOO was detected 
via the SZ effect in the first 10 months of the Planck survey and 
confirmed by XMM-Newton X-ray observations ( [Planck e'er 
laboration IX[[2011[ ). The cluster redshift was estimated by op 


tical griz observations performed with the 0.80-m IAC80 tele¬ 
scope, and through X-ray spectroscopic observations of the Fe 
K line, gi ving Zphot = 0.38 + 0.04 and z pe = 0.31 + 0.04 re¬ 
spectively ( [Planck Collaboration IX|201~T] ). The cluster mass es¬ 
timate, derived from XMM-Newton X-ray observations via it¬ 
eration about the M 500 - Yx relation of[Arnaud et al.[([2010|), is 
M 500 = (5.60 + 0.22) X 10^4 and the [0.1-2.4] keV X-ray 
luminosity was L 500 = (3.36 + 0.08) x 10"^^ erg s“^. 

Here we use deep (Ihr) r-band and shallower (~ 30 min) 
giz images obtained with SuprimeCam at the Subaru telescope, 
to derive the mass of the cluster by weak lensing. The paper is 
organized as follows. A summary of the observations and data 
reduction strategy is presented in Sect. Techniques adopted 
for star/galaxy separation are described in Sect.j^ A comparison 
with cluster properties derived by optical photometry is given 
in Sect. The weak lensing analysis is discussed in Sect. 
We attempt to refine the hydrostatic mass estimate from XMM- 
Newton data in Sect. Summary and conclusions are given in 
SectCl 

A standard concordance cosmology was adopted: = 0.7, 

Qm = 0.3, Hq = 70 km s“^ Mpc“\ giving a scale of 4.7 
kpc/arcsec at the redshift of PLCK GlOO.2-30.4 (z = 0.36, as 
derived in Sect. 0 - 


2 . Observations and data reduction 

PLCKGIOO (Fig. □ Fig. 13 was observed in the second part of 
the night on July 23, 2012 in the griz bands, with the wide-field 
Suprime-Cam camera ( [Miyazaki et al.|20(}2| ) mounted at the 8 m 
Subaru telescope. The camera is composed of 10 CCDs, each 
with 2048x4096 pixels and pixel scale 07207, allowing coverage 
of a field of 34'x27'. A sequence of 180s exposures was obtained 
following the dithering pattern recommended for SuprimeCam. 
The deepest (Ihr) image was taken in the r-band, to be used for 
the lensing analysis, with optimal seeing conditions (~ 075). Im¬ 
ages in g and i were obtained with a total exposure time of 1800s 
as planned, while due to technical reasons only 8 out of the 10 
planned exposures could be done in z, for a total exposure of 
1440s. Details of the observations are reported in Tabj^ 

The prereduction (overscan correction, bias, flat-fielding and 
masking of bad columns and autoguider) wa s based on the SD- 
FRED2 ( [Yagi et al.|2002[[0uchi et al.|2004[ ) pipeline developed 

for Suprime-Cam. _ 

The tool AstromC ( [Romano et al.[| 20 T 0 | ) was then used to 
compute the astrometric solution describing the deformations in 
the field of view of individual exposures. For each filter, the as¬ 
trometric solution was computed taking the NOMAD catalog as 
reference, and at the same time minimizing the differences in the 
position of the same sources measured in different exposures: 
this allowed us to obtain an internal accuracy of ~ 0.01 arcsec. 
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Table 1. Summary of the observations. The average seeing measured on the coadded images, and limiting magnitudes for pointlike sources at 
SNR=5,10 are also reported. 


Date 

Band 

Number 

Total time 

Average seeing (arcsec) 

mag((T =10) 

mag(o- = 5) 

23/07/2012 

g 

10 

1800s 

0.7 

24.9 

25.6 

23/07/2012 

r 

20 

3600s 

0.5 

25.5 

26.5 

23/07/2012 

i 

10 

1800s 

0.6 

24.7 

25.2 

23/07/2012 

z 

8 

1440s 

0.6 

23.8 

24.7 


better than the accuracy of the NOMAD catalog (rms ~ 0.2 arc- 
sec). Images wer e then resamp led and finally coadded using the 
SwARP softwar^ ( Bertin|2Q10 ). 


Photometric calibration was done by observing a field se¬ 
lected from the Stripe 82 area, and using the SDSS Stripe 82 
Standard Star Catalog ( [Ivezic et al.||2007| ) to derive photomet¬ 
ric zero points and SuprimeCam vs. SDSS color terms. Lim¬ 
iting magnitudes for point-like sources at signal to noise lev¬ 
els of SNR=5,10 were computed as in [Radovich et al.| ( |2004| ): 
simulated images with the same depth and background rms as 
each science image were produced using SkyMaker|^ software, 
and stars added at random positions; errors on magnitudes were 
computed from the difference between input and measured mag¬ 
nitudes, and used to derive the median magnitude for each value 
of SNR. 

In addition to the stacked image obtained as described above, 
we also produced a set of stacks where the PSF was homoge¬ 
nized to the same gaussian shape for all filters. This was done as 
follows: 


1. For each exposure, all CCD images were first regridded and 
combined in one image using SWarp, adopting the same cen¬ 
ter, pixel scale and size as the coadded r-band stack derived 
above. 

2. For each of them, stars were identified and the PSF fitted as 
a function of pixel coordinates using the PSFex softwar^ 
( |Bertin|2()TT] ), as described later in Sect.j^ 

3. The PSFMATCH task in IRAF was used to compute the ker¬ 
nel that transforms that PSF into a Gaussian with the tar¬ 
get FWHM, chosen to be 0.9 arcsec, the worst seeing value 
measured in all exposures and bands: since the PSF changes 
across the field, we divided each image in sections of size 
100 X 100 pixels, and computed the kernel in each of them. 

4. For each filter, all convolved exposures were finally summed 
together. 


3. Star-galaxy classification 

The selection of non-saturated stars and galaxies was under¬ 
taken in the magnitude (MAG_AUTO) vs. size (S = MU_MAX- 
MAG_AUTO) space ( [Huang et al |2QTT] ). MU_MAX is the peak 
surface brightness above background, which is constant for sat¬ 
urated stars. As in [Huang et al.[ ( [20lT] ), we rejected those sources 
with S lower than for stars as spurious detections. 

As a test of both the photometric accuracies and of the star- 
galaxy classification, we compared the g - r and i - z observed 
colors of stars with those derived at the Galactic plane position of 
PLCKGIOO by the TRILEGAL code ( [Girardi et al.|2Q05l[MT2l ). 
TRILEG AL allows the simulation of broad-band photometry of 
stars in any Galaxy field. The effect of the Galaxy dust extinc¬ 
tion on stellar colors was also included in the models, starting 
from the average and rms E(B - V) values derived from the 
Schlegel maps, and assigning to each simulated star a random 
value within this distribution. The transmission curves obtained 
by the combination of SUBARU/SuprimeCam CCD, telescope, 
atmosphere and filter contributions were adopted. Model and ob¬ 
served stellar colors are displayed in Fig[^ The zero point offsets 
required to optimize the overlap of observed and model colors 
are: -f0.038 mag (g), -0.030 mag (r), -0.023 mag (/), -0.016 mag 
(z). 

Stars with SNR> 100 were extracted and used in the PSFex 
software, designed to model the PSF variations in astronomical 
images. The code allows to model by a polynomial the pixel-by- 
pixel PSF intensity in a stamp of a given size as a function of 
the position in the field of view. Starting from this model, we 
can build a PSF stamp at each galaxy position, and derive the 
terms needed for the PSF correction of ellipticities, as described 


in Sect. 5.1 Different combinations of polynomial orders and 
stamp sizes were tried; based on the analysis of the residuals 
of the correction of the anisotropic PSF component (see later in 
Sect. in), we selected a polynomial of order 6 and stamps with 
size of 20x20 pixels (~ 5 x FWHM). 


These PSF-homogenized stacks were used for the photomet¬ 
ric measurements only, not for shape measurements, for which 
we used the unconvolved r-band image. A multi-ba nd photo- 
metric catalog was extracted by running SExtractoe[^ ( Bertin &] 


Amouts|1996[ ) in dual mode, taking the r-band as the detection 


image, and measuring the photometry on the convolved stacks. 

The Galactic dust extinction in the field cove red by Suprime- 
Cam was computed using the Schlegel maps ( [Schlegel et al.[ 
1998[ ), obtaining an average value and rms: £’b-v = 0.12 + 0.01. 


For each galaxy, observed magnitudes were corrected in each 
filter by adopting the dust extinction computed at its position. 


4. Cluster photometric properties 

Figure displays the composite gri image of the cluster field: 
the brightest 10 galaxies within 1 arcmin from the X-ray peak 
( [Planck Collaboration IX[[20TT] ) are marked. Two galaxies of 
similar brightness (r ~ 18.20 mag) are present in this region. 
We select as the center of the cluster the galaxy that lies closest 
to the X-ray peak, and identify it as the Brightest Cluster Galaxy 
(BCG). Several strong-lensing gravitational arcs are also visible, 
the brightest of which lies at ~ 40 arcsec from the BCG. Another 
lies at ~ 60 arcsec from the BCG. 


^ http://www.astromatic.net/software/swarp 
^ http ://www. astromatic.net/software/skymaker 
^ http ://www.astromatic.net/software/psfex 
http ://www. astromatic.net/software/sextractor 


4.1. Red sequence 

We selected the galaxies located in the inner 5 arcmin around 
the BCG that were classified as cluster members according to 
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Fig. 3. Stellar colors {top: / - z vs. g - r; bottom: r - i vs. g - r): mea¬ 
sured data (red dots) are compared to the values obtained by TRILE- 
GAL (green dots). 
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Fig. 4. rvs.g-r color plot of galaxies in the field. Galaxies classified 
as likely cluster members from their colors are plotted as green dots; 
the red dot is the BCG galaxy. The best fitting red sequence line is over¬ 
plotted. 


their colors, as described later in Sect. |5.2| The red sequence, 
g-r = a-\-br, was fitted by a biweight algorithm, giving a = 2.6, 
b = -0.05 (Fig.ig. Early-type galaxies were selected as those 
lying within 2cr from the best fitting line. Photometric redshifts 
were computed using the Zebra code ( [Feldmann et al.|[2006| ): 
the cluster redshift, derived from the average photo-z of r < 22 
mag cluster members, isz = 0.36 + 0.08. For comparison, we 
obtain z = 0.41 ±0.10 from the Sloan Digital Sky Survey (SDSS) 
photometric redshifts (DR9). 


4.2. Richness 


N 200 , the richness of the cluster within r 2 oo, was derived both 
from the selection of red-sequence cluster galaxies, and via mak¬ 
ing a statistical background subtraction. We used the r 2 oo de¬ 
rived from the weak lensing analysis detailed below in Sect. 
r 2 oo ~ 6 arcmin. Absolute r-band magnitudes and k-corrections 
were computed for all galaxies with the cluster redshift, z = 0.36, 
and to a magnitude limit of Mabs = -20 mag. The richness derived 
from selection of red sequence galaxies produces N 200 = 150. 
The richness estimate based on statistical background subtrac¬ 
tion was done following |Popesso et al.| ( |2005| ): the density of 
background and foreground galaxies was then computed in each 
of 20 sectors in an annulus with thickness 0.1 deg, starting at 
^200+0-1 deg. Sectors with too few galaxies were rejected, and 
the average density was finally computed, giving the number of 
background/foreground galaxies in the r 2 oo area of the cluster. In 
this way, we derived a richness N 200 -135. 


The richne ss can be translated to a mass estimate using 
e.g. Eq. 16 in |Andreon & Hurn| ( |2010| ): \ 0 gM 2 m = (0.96 + 
0.15)(logA^20o - 1.5) + 14.36 + 0.04. From these two consistent 
values of the richness, a mass \ 0 gM 2 m - (15.0 + 0.1) M© is 
derived. 
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5. Weak lensing analysis 

5.1. Shape measurement 


The derivation of the shear, and the correction of the effect of the 
PSF on galaxy shapes, is based on the Kaiser Squires and Broad- 
hurst approach ( [Kaiser et al.||1995[ KSB hereafter). In the last 
decade there has been a large effort to develop ne w algorithms to 
improve the accuracy in shape measurement (see jKitching et al. 
|2012| for a comprehensive review), which is crucial in particu¬ 
lar to derive cosmological parameters from weak gravitational 
lensing (cosmic shear). However, KSB is still widely used in the 
case of the analysis of massive structures such as galaxy clusters 
(see e.g.jOkabe et ST 2010t|Hoekstra et al.|2012]|Umetsu et al. 


2014t [Applegate et al.||2014|). O ur implementation follows that 

described in [Huang et al.[ ( [2QTT] ), with some differences that are 
outlined below. 

The KSB algorithm is based on the measurement of the 
following quantities for each galaxy, derived from the source 
brightness moments: the source ellipticity vectoij^^obs, and the 
smear polarizability {P^"^) and shear polarizability {P^^) tensors. 

The source ellipticity corrected for the anisotropic compo¬ 
nent of the PSF, ^aniso , IS first derived as: 


^aniso ~ ^obs 


- 


( 1 ) 


P = (2) 

where starred terms indicate that they are derived from measure¬ 
ment of stars. This is then related to the unknown intrinsic ellip¬ 
ticity e of the source and to the reduced shear, g = y/(l - K),hy\ 

^ In the following, the vector components along the X and Y axes are 
denoted as e.g. el and e2, respectively. 
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Saniso = e + P^g, 


(3) 


where the effect of the (isotropic) seeing is corrected by the term 
, dubbed by Luppino & Kaiser] ( |1 997 1 ) as the pre-seeing shear 
polarizability: 


psh* 

p7 psh psm_ _ psh t>sm ^ 


■ P^^q. 


(4) 


The tangential and cross components of the reduced shear, gt 
and gx are defined as: 

gt = -g\ cos lif - g 2 sin 2(f; g^ = -gi sin 2(p + g 2 cos 2(f, (5) 

<p being the position angle of the galaxy with respect to the as¬ 
sumed cluster center. 

The source brightness moments are weighted by a window 
function to suppress t he outer parts of th e galaxies which are 
noise dominated. As in |Huang et al.| ( |2QlT] ), the window function 
size was set to the value that maximizes the ellipticity signal to 
noise ratio SNe, defined as: 


SNe(6>) 


f I(0)W(l0l)d^0 

^ sky ^ f w^{\0\)cPe 


( 6 ) 


Sources with SNe < 5 were rejected from the analysis. We fur¬ 
ther excluded those galaxies with Py < 0.1, for which the PSF 
correction of the ellipticity is less reliable. 

The final output of the pipeline is the vector ^iso = ^aniso/^^- 
Provided that the average intrinsic ellipticity vanishes, {e) = 0, 
the average reduced shear is {g) = (^iso)- it is therefore common 
to set g = ^iso. _ _ 

As first pointed out by |Hoekstra et al.| ( [l998| ), it is important 
to adopt the same size of the window function to compute the 
galaxy ellipticities and the stellar correction terms used for the 
PSF correction. In [Huang et al.| ( [MTT] ) we computed the stellar 
terms in bins with 0 varying between the minimum and maxi¬ 
mum values allowed for the size of the galaxies (2 and 10 pixels 
respectively), with a step of 0.5 pixels. For each bin, they were 
then fitted by a polynomial, as a function of pixel coordinates, so 
that they could be derived at the position of each galaxy. This ap¬ 
proach presents some disadvantages: (i) it is necessary to make 
a polynomial fit for each PSF correction term and for each win¬ 
dow function size; (ii) large window sizes increase the amount 
of noise in the measurement of moments, even if we are using 
bright stars, which may contribute to different biases as a func¬ 
tion of galaxy size. A more rece nt approach consists of mod¬ 
elling the PSF ( jBerge et al.|2Q12| ), enabling us to reproduce the 
PSF at each galaxy’s positi on: this is now feasible with accura¬ 
cies ( Kitching e^al. 2013 ) o-{e) < 10“^. To this end, here we 
used the PSFex software, which allows modelling of the PSF 
as a linear combination of basis vectors and fits their coefficients 
with a polynomial as a function of the pixel position. For each 
galaxy position, PSF correction terms are then derived using ex¬ 
actly the same window function adopted for that galaxy. A sim¬ 
ilar approach was adopted in [Gruen et aL] ( |2013| ), who showed 
that using PSFex to model the PSF provides the accuracy re¬ 
quired for weak lensing analysis of galaxy clusters. The accuracy 
of the PSF fitting was checked by deriving the anisotropy cor¬ 
rected ellipticities for stars, and verifying the absence of system¬ 
atic effects in the residuals across the field of view (see Fig.[^. 

^ http ://www.astromatic.net/software/psfex 


Measured 



Corrected 



Fitted 



2000 6000 
X 


e1a -6.426e-04 +/- 1.548e-04 
O e2a 1.644e-04 +/- 1.420e-04 



Fig. 5. PSF anisotropy correction: the first three panels show the ellip¬ 
ticity pattern (measured, derived from the PSFex model and corrected 
(^aniso). X and Y are in pixels). The scale is displayed by the arrows in 
the middle (^ = 0.1). In the last panel, blue and red dots are the values of 
^aniso derived on real stars and on the PSFex model respectively; green 
dots are the values after the correction. 


Using a polynomial of order 6 as described above, we obtained 
(^'laniso) = (-6 ± 2) X 10-^ <e2a„iso> = (2 ± 1) X Fig. § 

shows the tangential component of the average stellar elliptic¬ 
ity before and after the anisotropic correction (< 6 x 10“'^), at 
different distances from the cluster centre. All of the steps out¬ 
lined above are implemented in the lensR pipeline introduced in 
[Huang et~ar] ( |2Qll| ). 

For each galaxy , we adopted the same weight definition as in 
[Huang et al. [ ( 201 1") : 


where A^o ~ 0.3 is the typical intrinsic rms of galaxy elliptic¬ 
ities. The uncertainty on the corrected ellipticities, A^iso, was 
computed following a numeric approach. For each star used in 
the PSF fitting, we compute the difference between the correc¬ 
tion terms measured on the real stars and on the PSF model 
derived by PSFex, and fit them as a function of the position. 
This allows estimation of an uncertainty on the stellar correction 
terms at each galaxy’s position. N values of ^iso were then com¬ 
puted (e.g. A^=100), randomly drawing both the observed source 
ellipticity and the stellar correction terms from a Gaussian distri¬ 
bution built from their errors. The error on ^iso used to compute 
the weight was taken to be the rms of these values. 


5.2. The background sample 

The dilution bias introduced by unlensed galaxies in the sample 
used for the shape measurement can adversely affect the cluster 
mass estimate. To minimise the contamination of this sample by 


unlensed sources, we followed the approach described in Formi- 


cola et al.] ( [2014[ ). Using the galaxy catalog with photometry in 
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Fig. 6. Average tangential PSF ellipticity at different distances from 
the cluster centre, before (black triangles) and after (red circles) the 
anisotropy correction. The shaded area is plotted at = +5 x 10""^. 

griz bands, we identified background, cluster members and fore¬ 
ground populations in a color-color diagram derived by choosing 
the appropriate combination of colors to be investigated. 

Our selection in the color-color (CC) space was based on the 
simultaneous analysis of galaxy colors and of the shear signal of 
the galaxies selected as background sources. We used the photo- 
metric data from the Cosmic Evolution Survey (COSMOS ; |Ilbert| 
et al.|2009] ) to train the color selection. The accurate photomet¬ 
ric redshifts of the sample allow a careful characterization of the 
distribution of the galaxies belonging to different redshift ranges 
in CC space. At the same time, we tuned the color cuts to max¬ 
imize the amplitude of the tangential shear component (Eq. 0 
measured from the background sources. 

Starting from the i-z vs g-rCC plot of the COSMOS galax¬ 
ies, we defined criteria based on color and magnitude cuts allow¬ 
ing to identify each galaxy population in the PLCKGIOO images. 
To this end, we computed the distance d of each galaxy in the i-z 
YS g-r diagram from the line repr esented in Eig.[7]( upper panel). 
This line (derived as described in |Eormicola et al.| ( |2014| )) sepa¬ 
rates the over-dense region mainly populated by sources at low 
redshift, from the over-dense region mainly populated by sources 
at redshift higher than that of the cluster. 

The background population was selected according to the 
following criteria: 

- galaxies with: 24.5 < r < 26.5 mag; 

- or. d < 0, 21.8 < r < 24.5 mag, i- z> 0.1; 

- or. d> 0, 21.5 <r < 24.5 mag, i- z < -0.3. 

The foreground and the cluster member population were 
identified according to these criteria: 

“ foreground galaxies: J>0, 21.0<r< 24.5 mag, -0.1 < 
i - z < 0.1; 

- cluster member galaxies: J > 0, r < 23.5 mag, 0.65 < g-r < 
2.0, 0.1 < / - z < 0.5. 



-0.5 0.0 0.5 1.0 1.5 

I-z 



Fig. 7. Selection in CC space. Upper panel: CC diagram for PECK 
GlOO galaxies. Cluster members are displayed with orange dots and 
foreground and background galaxies with green and magenta dots, re¬ 
spectively. Unclassified sources are plotted with black dots. Bottom 
panel: radial number density profiles of galaxies. The background den¬ 
sity profile (magenta squares) shows a decrease in the central region; 
overplotted is the magnification effect expected for a NEW profile. 


The total sample of background galaxies is represented with 
magenta dots in Fig.(upper panel), while foreground and clus¬ 
ter galaxies are marked with green and orange dots, respectively. 
Radial number density profiles for each identified population are 
shown in Fig. [^(bottom panel). The radial number density pro¬ 
file for cluster member sources (orange triangles) indicates a 
clear identification of the cluster member population as shown 
by the high clustering at small radii. As expected, the foreground 
population (green diamonds) appears uniformly distributed. The 
background sample (magenta squares) shows a drop in number 
counts towards the cluster centre, which is expected due to the 
so-called magnification bias ( [Broadhurst et al.|2008| ). Also dis¬ 
played in Fig. 1^ is the depletion in counts expected for a NFW 
profile with the parameters derived by the weak lensing analysis 
in Sect. |5.4| showing good consistency with the observations. 
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Fig. 8. r-band luminosity-weighted map in the PLCKGIOO field. The 
overplotted contour lines are the weak lensing S-map showing the mass 
distribution derived by weak lensing, at levels cr = 3, 4.5, 6, 8.5. The 
cross marks the position of the BCG assumed as the centre of the halo. 
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Fig. 10. Wavelet reconstructed X-ray image in the [0.5-2.5] keV energy 
band, obtained from the combination of the three XMM EPIC camera 
images. The contours show the weak lensing S-map (see Fig[^, at lev¬ 
els cr = 3.5,5,7,9 levels. The circle is the cluster centroid derived by 
the XMM analysis. The cross and square mark the position of the two 
brightest cluster galaxies (1 and 2 respectively in the text). 


5.3. Simulations 


r (kpc) 

1000 2000 3000 4000 



0 (arcsec) 

Fig. 9. Shear profile detected for PLCKGIOO (tangential and cross com¬ 
ponents, where the latter is expected to be null), and for comparison the 
values derived adopting the best-fit models. We indicate with mnfw the 
model obtained using the [Duffy et al.| ( [2008] ) relation and with nfw the 
model obtained keeping both M^tr and Cv/r as free parameters. 


Using the photo-z derived by jllbert et ( |2009| ), we ex¬ 
tracted the redshift distribution of the COSMOS galaxies with 
color and magnitude cuts describing the background population, 
needed for the evaluation of the critical density in the weak lens¬ 
ing analysis, as described in Sect. |5.4| 


The ellipticities derived by most of the available shape m easure- 
ment algorithms are affected (see jHeymans et al.|2006| and ref¬ 
erences therein) by a multiplicative (m) and an additive (c) bias 
defined as: ^obs “ ^tme = ^ ^tme + o. The value of m depends 
on the PSF type, as well as on the brightness and size of the 
galaxies, and in the case of KSB it is typically around 10%-20% 
( jHeymans et al.|2006] ). 


To investig ate this, we produ ced simulations using the Gal- 
Sinj^ software ( [Rowe et al.|2014) , that was specifically designed 
to simulate sheared galaxies within the GREAT effort (see [Man-] 
jdelbaum et al.|2013| and references therein). GalSim offers the 
possibility to use both analytic models for the galaxies, and/or 
a sample of 26 0000 real galaxies from Hubble Space Telescope 
COSMOS data. In addition, it can take as input the model pro¬ 
duced by PSFex, and hence can produce an image with the same 
PSF shape and spatial variations as the real image. We built 
our simulations so that ellipticities were measured on individ¬ 
ual stamps of sheared galaxies, in order to avoid the effect of 
blending from neighbour sources. 

To evaluate the bias, we first produced images with constant 
shear values, g = 0.01,0.05,0.1,0.15,0.2, for each ellipticity 
component. We randomly drew 10000 galaxies from the full 
sample, and repeated the simulation 10 times. For each shear 
value, we then verified that the expected input value could be 
recovered and evaluated the bias as a function of signal to noise 
ratio, as well as at different positions in the image to check for 
systematics related to a poor PSF correction. From these tests, 
we derived a multiplicative bias m = 10%, and no significant ad¬ 
ditive bias; we were not able to detect a significant dependence 
on the signal to noise ratio. 


^ https://github.com/GalSim-developers/ 
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Table 2. Best-fit NFW parameters, at different overdensities, derived taking the two brightest galaxies as cent re (n = 1,2): for e ach centre, in the 
upper row (a) both Mvir and Cyir were taken as free parameters; in the next ro w (b) the (Mvir, Cyir, Zd) relation ( [Duffy et al.|2008| ) was used; in the 
last three rows we use the c - M - z 2D relations in Meneghetti et al. (2014 i (c: all, d: relaxed, e: super-relaxed). Masses are in units of IO^'^Mq, 
radii in kpc. 
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Secondly, we used GalSim to simulate the shear signal pro¬ 
duced by a cluster having a NFW profile with redshift, virial 
mass (M ~ 10^^ h“^ Mq) and concentration parameter (c 2 oo = 4) 
similar to what we find here for PLCK-GIOO. We made ten 
such simulations and fitted the results in the same way as in 
Sect. |5.4| where the only free parameter was the mass and the 
concentration parameter was set to the input value of the simu¬ 
lation. The mass obtained applying the multiplicative bias was [Brainerd 
M 200 = (1.0 ± 0.2) X h-^Mo. 


mass estimate. As described in Sect. |5.2| the value of (fi) was 
computed using the photometric redshifts from the COSMOS 
sample, applying the color and magnitude cuts adopted for the 
selection of background galaxies. We find (J3) = 0.54, equivalent 
to assuming that all galaxies lie at the same redshift = 1.04. 

The shear profile was first fitted assuming a NF W halo and 
using the analytical expression for the shear given by [Wright &| 


([2000[): a Maximum Likelihood fitting approach was 


adopted ( [Radovich et al. [[2008] [Huang et al.[[2011[ ). In order to 


5.4. Mass estimate 

For the weak lensing analysis we used the sources in the cat¬ 
alog satisfying the following criteria: (i) classified as extended 
(Sect. 1^; (ii) with Py > 0.1, SNRolO; eliso, e2iso < 1 ; (iii) 
outside regions masked due to bad columns, haloes or spikes 
nearby bright, saturated stars. We further considered only those 
sources classified as background galaxies in Sect. [5.2[ This pro¬ 
duced a catalog of 17 000 sources in a 0.227 sq. deg area, corre¬ 
sponding to a density of 21 gals arcmin"^. 

We first extracted (Fig.[^ the S-map (see [Huang et al.|2()TT 


and references therein), showing the spatial distribution of the 
weak lensing signal. The peak is consistent with the position of 
the galaxy that we identified as the BCG of the cluster in Sect.|4 
We therefore took this position as the center of the cluster, anc 
used it to compute the cross and tangential components of the 
ellipticities whose profiles are displayed in Fig.[^ 

Computing the mass from the shear requires an estimate of 
the critical surface density: 


cH47rGDi/3r 


where p = DiJD/, Dis, Dg, and D/ are the angular dis¬ 
tances between lens and source, observer and source, and ob¬ 
server and lens, respectively. Since we have no spectroscopic 
training set, we are not able to judge the reliability of our pho¬ 
tometric redshifts for the lensed background galaxies. We there¬ 
fore adopted the approximation of a constant p for all galaxies. 
This introduces an overestimate of the average shear by a factor 
l-^((y^^) / - i)a:~ l-r0.2A:atz = 0.36, as discussed in 


Hoek- 


[stra et al.] ( [2000[ ). The correction of this bias is included m our 
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see how the mass estimate is affected by the assumption on the 
cluster centre and on the model fitting parameters, we performed 
different tests: {a) selecting both of the two brightest galaxies as 
the cluster centre (n = 1,2), and {h) trying both a 2-parameter 
fit (virial mass Myk and concentration Cyk), and by assuming a 
relation between Myk, Cyk and th e cluster redshift Zd - In the latter 
case, we first used the relation in [Duffy et al. ( 2008[ ): 




(9) 


with = 2.x /h Mq, A = 7.85, B = -0.081, C = -0.71. 


( 8 ) 


Then, we ado pted the c - M-z relations discussed in [Meneghetti 
[et al.[ P014[ see their Table 2). These were based on numeri¬ 
cal simulations of halos at different redshifts and halo relaxation 
states, undertaken for the analysis of the CLASH cluster sample. 
We used their 2D relations derived for NFW models for all ha¬ 
los, relaxed halos and super-relaxed halos. The results are given 
in Table The differences in the fitted parameters are within the 
uncertainties, with a slightly higher mass derived for ^ = 1. 

A model-independent estimate of the mass was then derived 
using the aperture densitometry method (see [Huang et al.|2011 


and references therein). This method provides an estimate of the 
projected 2D mass within a radius M 2 d(< ^m)- As discussed 
e.g., in [Qkabe et al.[ ( [20T0| ), before comparing the aperture (2D) 
and the NFW (3D) mass we first need to evaluate the correc¬ 
tion factor f(6m) = M 2 d(< 0m)IM2r>(< 0^) > 1. To this end, we 
used a gain the analytical expression given by [Wright & Brain^ 


( [2000 ) to compute fiPm) for an NFW halo with the best-fit val 
ues of Cyk, Myk found above (first row in Tab. [^. We obtained 
f = 1.39 at r 5 oo and f = 1.31 at ryk. The deprojected mass val- 
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Fig. 11. The radial profiles of the X-ray properties are displayed as follows. Top, left: surface brightness in the [0.5-2.5] keV range. Points and 
the solid line represent the measured data and the best-fit modified beta profile, respectively. Top, right: electron density. Solid line and grey 
area represent the best fit model and the relative error, respectively. Bottom, left: measured temperature profile. Points and solid line represent 
measured data and the temperature profile, respectively. Bottom, right: 3D temperature profile. Solid line and grey area represent best fit and the 
error, respectively. 


ues are: = (6 ± 2) x lO'^ Mq and = (1.0 ± 0.4) x lO'^ 

Me. 

The mass derived by weak lensing can be affected by the 
presence of large-scale structures along the line of sight ( |Hoek-| 
stra||2003j |Hoekstra et al.||2011| ). According to Figs. 6 and 7 in 
Hoekstra (2003), this uncertainty is of the order of ~2xl0i4ri 
Mq for a cluster with M = lO^^M© at z ~ 0.4, = 15 arcmin, 

which is comparable to the uncertainties derived in the fitting 
listed in Tabled 


6. X-ray analysis 

PLCK GIOO.230.4 has been observed for 13.7 ks with the 
XMM-Newton (obs-ID 0656202301) European Photon Imag¬ 
ing Cameras (EPIC), during the follow-up campaign of Planck 
newly detected clusters. This dataset was re-processed using the 
XMM-Newton Science Analysis System (SAS) version 14.0.0 in 
order to select the relevant calibration files. To remove any flare 
contamination, i.e. observation periods with an unusual count 
rate, we applied a 3cr clipping and a temporal wavelet Alter to 


the histogram of the photon arrival times, with a bin size of 50 s, 
in the energy band [1-5] and [10 - 12] keV. After the flare re¬ 
moval procedure we And a useful exposure time of 13.7 ks and 
5.2 ks for the two EPIC-MOS and the EPIC-PN datasets, respec¬ 
tively. To prepare the following imaging and spectroscopy anal¬ 
yses, we rebinned the flare-cleaned list of photon events in sky 
coordinates and energy, and associated both a 3D effective expo- 
su re and a background noise m odel to this data cube (see details 
in |Bourdin & Mazzotta|2008| ). 


The background noise was modelled after subtraction of all 
detected point-like sources using two components: the instru¬ 
mental noise caused by high energy particles interacting with 
the detectors, and several as trophy sical components, namely the 
Cosmic X-ray background (CXB) (e.g. [Giacconi et al.||2001|) 


plus a diffuse emission of Galactic origin (e.g. Snowden et al. 
|1995| ). For the instrumental component we used an analytical 
model inspired by the works of Kuntz & Snowden] ( [2008| ) and 
[Leccardi & Molendi| ( |2008| ), which holds for each EPIC camera 
a continuum plus several fluorescence lines (see detai ls in|Bour- 


din et al.|2013| ). As proposed by [Kuntz & SnowdEi| ( |2006| ), the 
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Fig. 12. Radial mass profile derived from the X-ray data, where the hy¬ 
drostatic mass profile and uncertainties are displayed by the solid line 
and the grey area. Error bars represent the 1 cr uncertainties described 
in the text. Also displayed are the mass profiles derived from the WL 
analysis (dashed line: NEW fit; points with error bars: aperture den¬ 
sitometry after deprojection). The inset shows M 2500 , ^ 1000 . ^500 and 
Mvir derived from our XMM analysis (blue squares); the Planck anal¬ 
ysis (green square); WL NEW, first row in Tab. (red triangles); WL, 
aperture densitometry (black circles): a small offset in radius is applied 
here for display purposes. 


Galactic foreground was modelled using two absorbed APEC 
( [Smith et~ar]|2001| ) spectra, with temperatures of 0.14keV and 
0.248 keV, respectively. The unresolved CXB spectrum follows 
a power law with L = 1.42 (Lumb et al. 2002| ). The particle back¬ 
ground was normalised in order to match the count rates of each 
camera outside the target region and in the hard energy band of 
[10 - 12] keV. The normalization of each astrophysical compo¬ 
nent was determined using a joint spectral fit outside any cluster 
emission, in the same region where the particle background was 
normalised. 


From the cleaned data we produced an image in the [0.5 - 
2.5] keV, which is background subtracted and corrected for any 
spatial variation of the effective area. From this image we deter¬ 
mined the X-ray peak at (RA = 350.56228, DEC = 28.52041). 
Fig. [T0| shows a wavelet denoising of this image. As we can see 
from this Figure, the X-ray peak (circle) lies within 5" of the 
galaxy that we assumed to be the BCG (cross), whereas the sec¬ 
ond brightest galaxy is 18" away. It is worth noting that there is 
also a good agreement between the direction of the elongation 
of the cluster, going from NE to the SW, as seen by the X-ray 
image and the isocontours from the lensing analysis. 

From the photon image we then extracted a point-source ex¬ 
cised surface brightness profile, centered on the cluster X-ray 
peak. In Fig. [m top left panel, we show this surface bright¬ 
ness profile, which we found to reach the background level at 
R = 1450 kpc. The temperature profile, shown in Fig. [TT] bot¬ 
tom left panel, was extracted in 4 bins from the same circu¬ 
lar region as the surface brightness The temperature was mea¬ 
sured in each bin by fitting the spectrum in the energy band 
[0.3-11] keV using an absorbed APEC thermal model with three 
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parameters, namely the normalization, temperature and metal- 
licity. The thermal models were convolved by the instrument re¬ 
sponse, as described in detail in jBourdin & Mazzotta|p008] ). The 
atomic hydrogen colum n density was fixed to its G alactic value 
Nh = 6.24 X 10^^ cm“^ ( Dickey & Lockman 1990| . 

To measure the X-ray mass we followed the "forward" pro¬ 
cedure described in [Meneghetti et al.j ( |2010| ) and [Rasia et ^ 


( 2012| ). In short, analytic profi les of the gas density and tem¬ 
perature jVikhlinin et al.j ( |2006|) were proje cted using the spec¬ 
troscopic like temperature [Mazzotta et al.j ( |2004| ) and fitted to 
the observed surface brightness and projected temperature. The 
3D information were then folded into the hydrostatic mass equa¬ 
tion (e.g. |Sarazin|[l988| or [Ettori et ar]|20I3| ). The confidence 
envelopes of the 3D profiles in Fig. |12| were computed fol¬ 
lowing a Monte-Carlo procedure. More specifically, we gener¬ 
ated several realisations of both a Poisson distributed surface 
brightness profile and a Gaussian distributed temperature pro¬ 
file, that hold the mean values and variances of the derived data 
quantities. For each of these realisations we fitted a new set of 
density and 3D temperature profiles, and stored these models 
when they obeyed the physical prior that the underlying hydro¬ 
static mass profile is increasing and positive. From = 500 
remaining hydrostatic density, 3D temperature and mass mod¬ 
els, we drew the confidence envelopes that encompass 68% of 
the realisations with minimal distance to their relative best fit. 
The envelope of the mass profile allowed us to derive a cluster 
mass of M 2500 = 3.2!0;« x Miooo = 4.4!°;^ x 

and M 500 = 5.4;!;^° X 10^"^M©, which corresponds to R2500 = 


541;|;^®kpc,/?iooo = Slh+^^kpc and/?5oo = 1093;!;“ kpc, re- 
spectively. The latter mass value is in excellent agreement with 
the mass derived from iteration about the Msnn - Yy relation , 
M 500 = (5.60 + 0.22) X 10^4 ( jpianck Collabomtion IX|2^ . 

Figure [T^ also shows a comparison of the X-ray total mass 
profile with those obtained from WL aperture mass densitometry 
and WL NFW model fitting (obtained with two free parameters). 
The inset shows the various mass estimates at density contrasts 
of 2500, 1000, and 500. The agreement between the different 
mass measurements is remarkably good at all density contrasts. 
Planck Collaboration XX| ( |2014| ) parameterised the differ¬ 


ence between their mass estimate and the true mass by the factor 
(l-b), where h contains all possible sources of bias, such as ne¬ 
glect of the non-thermal pressure support in the hydrostatic as¬ 
sumption, multi-temperature structure in the ICM, instrumental 
calibration uncertainties, etc. Assuming a weak lensing recon¬ 
structed mass of M 5 oo,wl ~ (7 + 1) x 10^"^ M© and comparing to 
the mass estimate obtained by [Planck Collaboration IX| ( |2011| ), 
we obtain (1-Z?) = 0.8 + 0.1 from the present observations. 

[Planck Collaboration XX[ ( [2014[ ) obtained a bias factor of 
{l-h) = 0.8 from their calibration of h using numerical simu¬ 
lations, and used this value, with a flat prior of [0.7,1.0], as 
the baseline for th eir cosmological analysis. [Planck Collabora- 


[tion XXIV[ ( [2015[ ) have updated this analysis to include a cali¬ 
bration of (1 - b) using weak lensing mass measurements, tak¬ 
ing into account sources of scatter such as triaxiality, the pres¬ 
ence of subtructures, and systematic errors in the weak lensing 
measurement itself. Assuming the weak lensing mass to be the 
true mass, their analysis yields (1 - /?)wtG = 0.69 + 0.07 and 


(1 - ^)cccp = 0.78 ± 0.09 for the samples of von der Linden 
et al.[ p014[ Weighing the Giants) and [Hoekstra et al.[ poT5 


Canadian Cluster Cosmology Project), respectively. Under the 
same assumption, the value we obtain from PLCKGIOO is com¬ 
patible within Icr with the bias factor derived from these updated 
measurements. 
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7. Summary and conclusions 


We presented here the first weak lensing analysis of the PLCK 
G100.2-30.4 cluster. The analysis confirmed that this is a mas¬ 
sive cluster (Myir > 10 ^^ Mq), in agreement with that expected 
from the richness estimate. 

At a density contrast of 500 with respect to the critical den¬ 
sity of the Universe at the cluster redshift, the weak lensing mass 
from aperture densitometry is = (6 + 2 ) x 10 ^^ M©, after 
deprojection. An NFW fit in the radial range 1-15 arcmin, shown 


by the dotted line in Fig. 12 yields a mass M 500 ~ (7 + 1) x 10 


14 


Mq. We checked the robustness of the weak lensing mass esti¬ 
mates by repeating the analysis with different choices of the clus¬ 
ter centre and of the NFW model fitting: in particular, we com¬ 
pared the results obtained by leaving both Cyk and Myir as free 
parameters, and using different (c, M, z) relations. These yield 
results consistent with those given above. 

Comparing the mass derived by weak lensing with the value 
obtained by the reanalysis of XMM data, we obtain a bias factor 
of (1-/?) = 0.8 + 0.1. This is compatible within Icr with the bias 
factor obtained by [Planck Collaboration XXIV| ( |2015| ) from their 
calibration of (I-/?) using newly-available weak lensing samples 
from |von der Linden et al.| ( |MT4| ) and |Hoekstra et al.| ( [2015| ). The 
present results illustrate the importance of high-quality data and 
the need for larger samples in such cluster mass studies. 
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